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SUMMARY 

Three semi-empirical aerodynamic stall models are compared with respect 
to their lift and moment hysteresis loop predication, limit cycle behavior 
prediction, easy implementation, and feasibility in developing the parameters 
required for stall flutter prediction of advanced turboprops. For the compar- 
ison of aeroelastic response prediction including stall, a typical section 
model and a plate structural model are considered. The response analysis 
Includes both plunging and pitching motions of the blades. 

In model A, a correction to the angle of attack Is applied when the angle 
of attack exceeds the static stall angle. In model B, a synthesis procedure 
Is used for angles of attack above static stall angles and the time history 
effects are accounted through the Wagner function. In both models the lift 
and moment coefficients for angles of attack below stall are obtained from 
tabular data for a given Mach number and angle of attack. In model C, referred 
to as the ONERA model, the lift and moment coefficients are given in the form 
of two differential equations, one for angles below stall and the other for 
angles above stall, lhe parameters of these equations are nonlinear functions 
of the angle of attack. The effects of vortex shedding, an Important feature 
of dynamic stall, are not considered in model A, accurately considered in model 
B, and approximately considered in model C. However, it is observed that the 
high frequency, low amplitude oscillations, sweep, and high subsonic Mach num- 
ber operating environment of advanced turboprops favor light stall conditions 
where the effect of vortex-shedding Is less severe. This permits the use of 
simple models like models A and C In the stall flutter analysis of advanced 
turboprops . 


INTRODUCTION 

Highly loaded propellers, called advanced turboprops, are proposed to 
power transport aircraft at high subsonic speeds. The renewed Interest In the 
propeller Is brought about by the significant benefits in fuel consumption. 
Flutter analysis of these turboprops is needed to determine the critical 
(flutter) speed below which the aircraft has to operate to avoid catastrophic 
failure. A response analysis is needed to determine the loads on the blades 
and fatigue life of a turboprop. Figure 1 shows a typical advanced turboprop 
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wind tunnel model, lhe blade is made thin to increase the drag divergence 
Mach number, and is swept to reduce the local effective Mach number. The sweep 
introduces a large degree of coupling between bending and torsion, which alters 
the aeroelastic behavior. In addition, to hold propeller diameter to a rea 
sonable value, and to have high disk loading, a large number of blades are 
provided for the advanced turboprop. Consequently the turboprop has to be 
analyzed accounting for aerodynamic interaction between the blades. The aero 
elastic response and flutter analysis of advanced turboprops is an ongoing 
research effort at NASA Lewis Research Center. 

Under normal conditions (ref. 1), the blade sections of a conventional 
propeller are at low angles of attack, the flutter speed is high, and generally 
there is an appreciable margin of safety between the operating speed and the 
flutter speed. However, during the take-off period, the propeller blade sec 
tions may operate at high angles of attack and could be subjected to stall 
flutter, a self sustained vibration in a periodically separated flow condition, 
lhe flutter speed at this condition is very low, and stall flutter affects the 
fatigue life of the turboprop. Hence, predicting this stall flutter speed is 
a critical design task, and an appropriate method of analysis of stall flutter 
is a critical research need. Prediction of flutter speeds at low angles of 
attack (classical flutter) for advanced turboprops has been performed by modi- 
fying the existing analysis methods as described in reference 2. Both ana 
lytical and experimental results are presented in reference 3. Recently, 
additional experimental flutter data was presented in reference 4 . lhe analy 
ses included the effect of number of blades (cascades) and the blade sweep. 
However, the studies on stall flutter are few and mostly empirical, lhis is 
due to the complexity in modeling the flow in a periodically separated flow 
state. 

lhe objective of the present effort is to develop stall flutter models 
for advanced turboprops. As a part of this general effort, the available 
dynamic stall models are reviewed and applied to simple structural models to 
study the extent of their validity, and to select an appropriate model for 
advanced turboprop application. In this report, three dynamic stall models 
are used together with typical section and plate structural models. Their 
performance is investigated from the view point of lift and moment hysteresis 
loop prediction, limit cycle prediction, easy implementation and the feasibil 
ity of developing Lhe corresponding stall models for arbitrary airfoils, lhe 
plate model in conjunction with its normal modes is considered in order to 
explore whether the analytical integration of the loads in the blade spanwise 
direction is feasible, lhe numerical study is performed for a single blade 
and cascade effects are not included. 
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static pitching moment curve slope at zero angle of attack 
semi-chord, m 

lift coefficient in linear range, equation (20a) 
lift coefficient in nonlinear range, equation (20b) 
lift coefficient given by linear relation, equation (20a) 
static lift coefficient 
unsteady c^, equation (18a) 

moment coefficient in linear range, equation (20d) 

moment coefficient in nonlinear range, equation (20e) 

moment coefficient given by linear relation, equation (20d) 

static moment coefficient 

unsteady c m , equation (19a) 

normal force coefficient 

chord length, m 

aerodynamic lift force coefficient 

aerodynamic moment coefficient about quarter chord 

aerodynamic moment coefficient about elastic axis 

empirical parameter, equation (19a) 

plunging degree of motion, positive downwards 

polar moment of inertia of airfoil mass about elastic axis 

torsional stiffness coefficient corresponding to pitching 
d i sp lacement 

bending stiffness coefficient corresponding to plunging 
d i splacement 

unsteady aerodynamic empirical factor, equation (3) 

reduced frequency, based on semi chord 

length of the plate, m 

Mach number, normal to leading edge 

mass of the airfoil per unit span 
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number of blade segments 
number of normal modes 

empirical parameters for C|_u, equation (18b) 
total aerodynamic moment about elastic axis, positive nose up 
total aerodynamic lifting force, positive downwards 
empirical parameters for C |_u , equations (18d) 

Reynold's number 

•^I a /mb2, radius of gyration about elastic axis in semi-chords 

empirical parameter, equation (20b); radial distance on the blade 

airfoil static moment about elastic axis 

nondimensional time measured from Instant of stall onset 

predicted value of s m when vortex leaves trailing edge, 
equation (16) 

reduced time, Vt/b; also empirical parameter, equation (20a) 

thickness of the airfoil; time 

wt, nondimensional time parameter 

airfoil thickness to chord ratio 

time when dynamic stall first occurs 

nondimensional flight speed parameter 

resultant velocity = Vn^r^ *■ v£, m/sec 

free stream velocity, m/sec 

S/mb, distance between the elastic axis and the mass center 
measured in semi chord, positive towards the trailing edge 

instantaneous angle of attack; pitching degree of motion, 
positive nose up 

amplitude of oscillation, deg 

dynamic moment stall angle, equation (15), deg 

effective angle of attack, deg, equation (1) 

decay parameter, equation (3) 
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static stall angle, deg 
Prandtl- Glauert number 

empirical constant, equation (18e) normally equals 0.18 

stall delay function, equation (2) 

functions defined in equations (4) and (5) 

shift in angle, equation (18b), deg 

shift in angle, equation (18c), deg 

Incremental dynamic stall angle, equation (1) 

difference between static lift and extended linear lift 

incremental lift coefficients, equations (18d) and (18e) 

difference between static moment and extended linear moment 

incremental moment coefficient, equation (19b) 

dynamic parameters, equation (18g) 
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empirical parameters, equations (16) to (18) 

empirical parameters for C,.,,, equation (19b) 

MU 

empirical parameter, equation (20a) 
m/irpb^ airfoil-air mass ratio 
nondimensional plunging displacement 
free stream air density 
damping parameter in plunging motion 
damping parameter in pitching motion 
empirical parameter, equation (20a) 

Wagner function, equation (14) 
rotation speed, rad/sec 
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frequency of harmonic oscillation, rad/sec 
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Stall, Dynamic Stall and Stall flutter 


As the angle of attack of an airfoil increases, the lift coefficient 
starts increasing. After a certain value of the angle of attack is reached, 
however, the lift drops suddenly because the flow over the airfoil separates. 
This is the condition of stall. Ihis separation process depends on the air 
foil shape (leading edge radius), thickness, Reynold's number, maximum thick- 
ness position, and Mach number. It has been recognized (refs. 6 and 6) that 
there are three principal types of stall: (I) trailing edge stall, where 

there is a gradual loss of lift at high lift coefficient as the boundary layer 
separation progresses gradually forward from the trailing edge, (2) leading 
edge stall, where there is an abrupt loss of lift, as the angle of attack for 
maximum lift is exceeded, with little or no rounding over of the lift curves; 
and is associated with the bursting of a laminar leading edge separation bub 
ble, and (3) thin airfoil stall, where there is a gradual loss of lift even at 
low lift coefficients, and develops when a separation bubble originates near 
the leading edge and lengthens progressively as the angle of attack increases, 
lhe behavior of the lift coefficient in each type of static stall is shown in 
figure 2. It has also been found that trailing edge stall occurs for airfoils 
having thickness to chord ratio t/c greater than 0.16, leading edge stall 
occurs for airfoils having t/c of 0.09 to 0.16, and thin airfoil stall for 
airfoils having t/c less than 0.09. 

lhe term dynamic stall refers to the unsteady separation and the stall 
phenomena of an airfoil oscillating into and out of stall, figure 3 (from 
ref. 7), shows the flow field structure, as well as the normal force and pitch 
ing moment characteristics throughout an oscillation cycle of an airfoil during 
dynamic stall for an NACA 0012 airfoil, lhe predominant feature of dynamic 
stall (ref. 7), is the shedding of a strong vortex like disturbance from the 
leading edge region, which alters the chord wise pressure distribution. This 
vortex moves down stream over the upper surface of the airfoil at about 36 to 
40 percent of free stream velocity, lhe unsteady aerodynamic forces due to 
the passage of this vorticity produce a lift and nose down moment, with values 
much greater than the corresponding static stall loads, lhe magnitude of the 
increase depends on the strength of the vortex and its distance from the sur 
face, lhe formation and movement of the vortex depends on the airfoil shape. 
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angle of attack and the rate of angle of attack. The distance of the vortex 
from the airfoil depends on the rate of angle of attack and the Instantaneous 
angle of attack. As the vortex leaves the trailing edge, a peak negative 
pitching moment is obtained. The airfoil then remains stalled until the angle 
of attack drops sufficiently for reattachment of the flow to occur. 

It can also be seen that airfoils during dynamic stall exhibit large hys 
teresis loops in both lift and pitching moment curves when viewed as a func- 
tion of angle of attack. This means that an airfoil with positive a stalls 
at an angle a greater than the section static stall angle, a ss , known as 
stall delay, while the stall recovery during negative <i occurs at an angle 
less than a ss . It can also be observed that the pitching moment coefficient 
shows loops representing contributions to negative damping (clock-wise loops) 
and positive damping (counter clock wise loops), and a net negative damping 
may lead to divergent oscillations. It is to be noted that for oscillations 
that occur wholly below or wholly within stall, there is always positive damp 
ing. Only for oscillations about a mean angle of attack near static stall can 
the net pitch damping become negative. 

In summary, dynamic stall begins at an angle of attack greater than the 
static stall angle, followed by the shedding of vorticity from the leading and 
trailing edges. As the airfoil oscillates in and out of stall, the dynamic 
forces and moments show hysteresis and can attain values that are far greater 
than their static counterparts. The dynamic stall phenomenon and its effects 
vary depending on the airfoil shape, reduced frequency, mean angle and ampli- 
tude of oscillation, Mach number, Reynold's number, type of airfoil motion, 
sweep, and three-dimensional flow effects. 

An important difference between stall flutter and classical flutter is in 
the character of flow. Stall flutter occurs with partial (light dynamic stall) 
or complete (deep dynamic stall) breakaway of the flow from the airfoil during 
at least part of every cycle of oscillation. This is in contrast to classical 
flutter, where the flow is attached to airfoil at all times. The essential 
feature of stall flutter is the nonlinear aerodynamic reaction to the motion 
of the airfoil. Ihis nonlinear nature allows, in principle, the prediction of 
the final equilibrium amplitude of vibration (a limit cycle). Ihis is differ 
en t from classical flutter where only the stability boundary is usually deter 
mined. In short, stall flutter refers to a self-excited and self- sustained 
vibration in a periodically separated flow condition. Stall flutter also 
differs from classical flutter in that the torsional and bending frequencies 
are not necessarily close together even though both modes contribute to the 
stall flutter. 

It has also been observed in stall flutter (ref. 8) that (1) there is a 
sharp drop in critical flutter speed, (2) the flutter frequency rises towards 
the torsional frequency, (3) the motion is predominantly torsion (single degree 
of freedom flutter), and (4) the stall flutter speed reaches a minimum and 
rises until the flow is completely stalled. 


DYNAMIC SI ALL MODELING 


lhe complexity in modeling stalling is due to the following reasons: 

(1) flow separation and turbulence effects during part of the cycle of oscil- 
lation; (2) a new variable, the mean angle of attack a 0 , is introduced into 
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the determination of aerodynamic derivatives and automatically requires con- 
sideration of the effects of Reynold's number and airfoil shape; (3) the aero- 
dynamic derivatives can no longer be simply superimposed (as in classical 
flutter) i.e., the aerodynamic results of a pitching motion cannot be separated 
from those of a simultaneous translatory motion; and (4) rate of angle of 
attack effects must be included in the analysis. 

Earlier work on unsteady stalled flow has been either wholly experimental 
or consisted of empirical and semi-empirical modifications to classical 
flutter theory. These earlier studies have been based on the hypothesis that 
the decrease in flutter speed may be due to a decrease in aerodynamic torsional 
damping. Halfman, et al., (ref. 8) reviewed these methods and presented exper- 
imental data for lift and pitching moment in pure pitch, and in pure plunge. 
Sisto (ref. 9) presented a nonlinear mechanics approach to the problem of stall 
flutter, and verified his predictions with experiments for both isolated and 
cascade of blades. Schnittger (ref. 10), used data of Halfman et al. to inves- 
tigate stall flutter in compressors and found that cascade effects have a sup- 
pressing effect on stall. Amer and La Forge (ref. 11) developed a procedure 
for calculating blade bending moments, torsion moments, negative aerodynamic 
damping, and the lift hysteresis during stall using Halfman, et al. experi- 
mental data. Baker (ref. 1) conducted experiments to measure the stall flut- 
ter speeds of thin wings representative of propeller blades. Rainey (ref. 12) 
measured aerodynamic damping to Investigate stall flutter. Ham (ref. 13) con- 
ducted experiments to investigate stall flutter of helicopter blades. Carta 
(ref. 14) used an energy principle to calculate aerodynamic damping, and used 
Halfman et al. data at high angle of incidence for turbojet engines. 

Recent research on dynamic stall followed two approaches, one theoretical 
(refs. 15 to 32), and the other based on experimental data, references 33 to 
41. These research efforts on dynamic stall are summarized in references 42 
to 47. The flow elements to be included in dynamic stall modeling are dis- 
cussed in reference 48. Reference 47 tabulated the dynamic stall prediction 
methods according to the technique used in the formulation, and roughly graded 
most of the models according to the salient features of each model. An updated 
version of the table is presented in table 1, from which the strengths and 
weaknesses of the models can be seen more readily. 

The theoretical approaches are the Navier Stokes methods, the discrete 
vortex methods, and coupled viscous inviscid methods (zonal methods). Navier 
Stokes methods (refs. 15 to 17) attempt to solve the relevant equations in 
their fundamental form by numerical techniques. The discrete vortex approach 
(refs. 18 to 24) normally ignores the viscous terms in the basic equations and 
assumes potential flow without the boundary layer. The viscous nature of the 
flow is modeled or taken into account, by the generation and subsequent induced 
transport of discrete combined vortices. The manner and location of their 
generation is normally obtained empirically or via appropriate boundary layer 
calculations. In the zonal methods (refs. 24 to 32) the various regions of 
flow, viscous, nonviscous, and transition regions, are modeled separately. In 
the numerical implementation of the model, the regions interact in an iterative 
manner. These theoretical models require a lot of computer time and are lim 
ited by the assumptions and restrictions of the formulation. So they are not 
suitable in a routine aeroelastic analysis. 

Semi - empi rica I methods (refs. 33 to 41) based on experiments attempt to 
simulate the gross features of stall. They have gained much interest for the 
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following reasons: (1) they attempt to use static data with a dynamic correc- 

tion. This is advantageous because static data can be easily generated and 
automatically includes the effects of Reynold's number, Mach number, and air- 
foil shape. (2) They take less computer time so they can be used in a routine 
aeroelastic analysis. 

In the published literature there are four types of models for dynamic 
stall based on experimental data ( semi- empi rical ) . The first (ref. 33) is to 
use measured data tabulated in a three-dimensional array (angle of attack, 
reduced pitch rate, reduced pitch acceleration) with an appropriate correction 
for local Mach number. This requires a large amount of data storage for each 
airfoil, frequency of oscillation and the associated interpolation. The second 
(refs. 34 and 35) is to utilize a corrected angle of attack when the angle 
exceeds the static stall angle. This correction is a function of the rate of 
change of angle of attack.* The third (refs. 37 and 38) is to reduce the 
large volume of data obtained from experiments to compact expressions (synthe- 
sization). The fourth (ref. 39) is to describe the lift and moment coeffi 
cients in terms of ordinary differential equations (ONERA model). Ihe last 
two models attempt to fit the experimental data by way of equations. It should 
be noted that in the first two models the time history of the motion is not 
taken into account. 

Ihe empirical parameters used in the semi - empi rical models referred above 
are usually obtained from airfoils oscillating in pitch about quarter- chord . 
However, an arbitrary motion includes both pitch and plunge motions. So stall 
modeling should include separate h and h terms in the Identification of 
the empirical parameters. This would seem to require experiments that allow 
for plunge as well as pitch motions. But, it would certainly be convenient if 
only dynamic-pitch experiments were needed and the plunge effects could be 
derived from them by other means. Many of the existing dynamic stall modeling 
methods do not take into account this distinction. Reference 49 maintains 
that the equivalence between pitch and plunge is not valid since the way the 
stall cell forms and propagates is different in pitch and plunge. However, 
tests have shown that the equivalence between pitch and plunge is reasonable 
in light stall but not in deep stall. For more discussion on the equivalence 
of pitch and plunge see references 50 and 51. In the semi empirical models, 
mentioned above, the ONERA model attempts to take this distinction into 
account. 

In the following the last three semi - empi rical models, which are desig- 
nated as A, B, and C, are described. These methods are selected on the basis 
of their availability and easy implementation. Equations are presented to 
calculate lift and pitching moment coefficients. Equations similar to those 
given for pitching moment coefficient can be used to calculate coefficient of 
drag. It is to be noted when calculating the resultant velocity on the air- 
foil that the present study is restricted to nonrotating structural models, 
hence the resultant velocity is same as the free stream velocity. 


*The Mil model of reference 36, corrects the angle of attack as a function of 
rate of angle of attack, but idealizes the loading due to vortex as an Impact load 
at the instant of stall occurrence. The magnitude of this Impact loading is 
obtained from experiments. 
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MODEL A - CORRECTED ANGLE OF ATTACK APPROACH 


This model is presented in references 34 and 35. In this model, the 

actual angle of attack, a, is corrected to obtain an effective angle of 

attack a r . Then, the coefficients c„ and c . . , which are functions of 
E 8. mc/4 

Mach number, are obtained from static airfoil data. The correction to the 

angle of attack is a function of rate of angle of attack and is based on 

oscillating airfoil test data. 

The relation between the actual and the effective angle of attack is 
given in terms of an Incremental dynamic stall angle, Aa DS , as 


Then an empirical relation is established from experiments between Aa^ and 
the rate of angle of attack as. 


Here V is the resultant velocity, c is the chord of the airfoil, <i is the 
rate of angle of attack, A is the nondimens ional rate of angle of attack, the 
function y has to be determined empirically from oscillating airfoil test 
data, and K-| is given by 


The function y is a function of Mach number, airfoil maximum thickness 
to chord ratio and has different values for lift and moment. This indicates 
that the Influence of y is not the same for both lift stall and moment stall. 

From equation (2), it can be seen that y is the slope of the straight 
line representing the variation of Aa^ s with -s/|A| . 1 he slope is constant 

for thicker airfoils. However, for thin airfoils, experimental data reveal 
that there are two slopes y^ and y^, depending on the values of ■>/ |A| . The 
value of ->y | A | at which the slope changes is designated as Break ' 

definitions of y^ , y^, Act^, and -v/ 1 A | are schematically shown in figure 4. 
The values of y, for both lift stall and moment stall are listed below: 


a £ = a - K i Aa DS 


( 1 ) 



( 2 ) 


K] = 3/4 + 1/4 sign(<i) 


(3) 


Lift Stall 


Y 2L = 1.4 - 6.0 (0.06 - t/c) 


If Mach No. < 0.4 + 5.0 (0.06 - t/c) 


If Mach No. > 0.9 <• 2.5 (0.06 - t/c) 


Y 1L - 0.5 


( 4 ) 
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Moment Stall 


Y 2m = 1.0 - 2.5 (.06 - t/c). If Mach No. < 0.2 


t 2M = °* 


If Mach No. > 0.7 + 2.5 (.06 - t/c) 


YlM = 0.0 


( 5 ) 


where subscripts L and M denote lift and moment respectively. It Is to be 
noted that In the above expressions, the effect of Mach number Is Included In 
the definition of and y^. The variation of y ^ with Mach number Is 

schematically shown In figure 4. Also based on experimental data, >/|A| 

Is expressed In terms of t/c as rea 

^Break = 0-06 + 1 -5(0.06 - t/c) (6) 


The next step Is to express Aa QS In terms of y^ , y ? , and \[\h\. 

%s ' <T U L.H) >4*1) (sign A) if 4m < 4m 6reak (7) 

A “0S ' ( t 1(L,M) '^^Break * T 2(I.,M) ( " '^Break , ) s ’ 9n (“) 

If 4m > ViAT Break (8) 


This formulation has the effect of reducing the angle of attack for positive 
d, and of Increasing the angle of attack for negative d. The final lift and 
moment coefficients are calculated as 


c 


8. 


\ at * L " % : 



(9) 


C ,. = C ( a_ „) 

mc/4 m v E,M 


( 10 ) 


where a Is the angle of attack for zero lift. 

s - 0 

It should be noted that in this formulation the effects of Mach number 
and thickness to chord ratio are explicitly Included. This will be useful in 
the application of this model directly to advanced turboprop blades, which are 
very thin, operate at different Mach numbers along the radius, and have a 
varying thickness to chord ratio along the radius. It should be noted that 
above formulation is used only when the actual angle of attack is above the 
static stall angle. This implies that after separation of the flow, the 
reattachment is assumed at static stall angle. Reference 38 observed that 
this may not reproduce the pitching moment loop correctly, thereby affecting 
the aerodynamic damping calculations. 
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MODEL B - SYNTHESIS PROCEDURE 


Reference 38 presents a synthesized approach to dynamic stall modeling. 

The synthesized data essentially consist of seml-empi rical ly obtained analyti- 
cal expressions representing qualitative approximations to the various observed 
physical features associated with the dynamic stall of airfoils. The model 
adequately accounts for the effects of formation and streamwise travel of the 
dynamic stall vortex. The process Involves curve-fitting analytic expressions 
to the test loop data with the objective of determining the unknown parameters 
or coefficients embedded in the analytical expressions. 

Three dynamic parameters are defined to predict the dynamic stall events: 
(1) the instantaneous angle of attack, a; (2) the nondimensional rate of angle 
of attack, A; and (3) a decay parameter, a w , which accounts for the time 
history effects of the change in a, and is based on the Wagner function, the 
data on which this synthesis procedure is based accounts for some effects of 
airfoil shape, Mach number, and Reynold's number. In the following, the decay 
parameter, a w is defined, then three stages of dynamic stall are identified 
and expressed in terms of the empirical parameters obtained from the synthesis 
procedure. The relations for unsteady lift and moment based on the above 
parameters are presented. 


Definition of a w 


For a two-dimensional airfoil going through an arbitrary change in angle 
of attack, one can describe an instantaneous effective angle of attack, a F , by 
using Duhamel's Integral (ref. 40), 


a E (s) 


a(0)4» c (S,M) 



- ».M) 


do 


( 11 ) 


where a(0) corresponds to the initial angle of attack, M represents Mach 
number, <$> c (s,M) (Wagner function) is the response to a step change in a, 
and s is the nondimensional time given by 

s = 2 V t/c (12) 

The decay parameter is given by 

a w = a(s) - a t (s) (13) 

The decay parameter, a w , is the difference between the instantaneous 
angle, a (s), and the effective angle , and, therefore, acounts for the 

time history effects of the change in a. The effect of compressibility is 
also Incorporated in the definition of a w , by defining 

. u , n ,,, - 0.0456 (1 M 2 ) _ -0.35 (1-M 2 )// 

<t>£ (s,M) = 1 . - 0.165 e - 0.335 e I 1 _ ^2 ^ (14) 

The calculation of A, and a w for sinusoidal motion and arbitrary 
motion are given in the latter sections. 
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Prediction of Dynamic Stall Events 

In the published literature, three stages of dynamic stall are defined: 
(1) the onset of stall, which gives the angle at which stall starts (see 
fig. 3); (2) the vortex at the trailing edge, which gives the time (nondimen- 
sional) taken for the vortex to travel to the trailing edge; and (3) the 
reattachment, the angle at which the flow reattaches to the airfoil. In the 
present model these three stages are related empirically by the nondimensional 
rate of angle of attack. A, the decay parameter, a w , and the static stall 
angle, a ss . The values of A, a w , and a ss correspond to the values at the 
instant of moment stall, since moment stall occurs prior to lift stall. 


Onset of stall: The angle of attack at which the stall begins is given by 

(16) 


a n = ( 1 . -t- e *• C A A + C a ) a 
Dm Am m wm wm ss 


where subscript m refers to values at 
empirical parameters e, C. , and C 
mental data. 


the point of moment stall and the 
are obtained from curve fitted experi- 


Vortex at the trailing edge: After the occurrence of moment stall, there 

is a significant increase in negative pitching moment due to the travel of the 
stall vortex, lhe maximum negative pitching moment occurs when the vortex is 
near the trailing edge of the airfoil, lhe time at which the vortex leaves 
the trailing edge is given by 


s . = 1.0/(C a ,A + C .«*_ ) 

mt At m at Dm 


(17) 


s . is the total nondimensional time for the vortex to travel from the 
mt 


where 

leading edge to the trailing edge, and the empirical parameters C 


At 


and 


C , are obtained from experiments 

at 


Reattachment: For Mach numbers < 0.4 , the reattachment occurs at an angle 

a ut which is less than the static stall angle. At higher Mach numbers, a 
K L K t 

can be greater than the static stall angle 


“Rl 3 (1 ' 1 

lhe empirical parameters C. ri and 

AR 

f rom experiments . 


*■ C A C a )a 
AR m wR wm ss 

C „ are curve fitted 
wR 


(18) 

for a given airfoil 


In summary, there are seven empirical parameters to predict the dynamic 
stall events and they depend on Mach number, Reynold's number, sweep angle and 
airfoil shape. 

Now the unsteady lift and moment coefficients, including dynamic stall 
effects, are expressed in terms of additional empirical parameters and are 
given below. 


Unsteady Lift Coefficient 


C l = C LU = C LS (a ' Aa l ’ Aa 2 } ' a ol Aa l + AC L1 + AC | 2 


(18a) 
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(18b) 

(18c) 


AC 


A “l = (P 1 A * P 2“w + P 3> a ss 
Aa 2 = A 2 a ss 

Li = Qi A * ^“w + ^3^ a/a ss^ * Q 4 (a/a ss )< 


(18d) 


AC 


L2 = Vl + Q 6 Aa 2 + Q 7 (a Dm ) 


1 - e 


OVm>' 


< B lV‘ 


(18e) 


2V <t - W 


m 


(18f) 


= 0 a < a 

“ SS 


-1 a < a < a n 

SS — Dm 


*i 



“Dm 1 


r, . M] 


_“ss ‘ 


\ s mt/ 


0 c s m < s . 
m - mt 


= 0 s > s . 

m mt 


0 

a. < 

a 

SS 



a _ 
a 

SS 

1 

“ss 

< a 

- a Dm 

*Dm 

a ss 

1 

0 < 

VI 

E 

on 

s mt 

E 

q! 

0 

- 1 

a - i 

“RE " 


a 

SS 

“TE 

“RE 

“re 


l 7E 


= 0 a < a 


RE 


(l8g) 


where 6] and depend on onset of stall and reattachment. 

In equation (18a), a 0 d is the conventional static lift curve slope. The 
shift in angle of attack, Aa-j and Ao^, are associated respectively with the 

unsteady effects below stall and with the occurrence of dynamic stall and 
reattachment. AC|j represents the unsteady effects over static lift for 
unstalled airfoils, and ACl 2 represents the effects associated with dynamic 
stall. Equation (18a) expresses the synthesized unsteady coefficient as a sum 
of static lift coefficient, C , at some shifted angle (a-Aa Aa ) plus an 
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increment lift coefficient (AC L1 + AC |2 ). For unstalled cases Ac^ and 

AC| 2 are zero. Also G-| = 0.18, s m is time measured from the instant of the 
occurrence of dynamic moment stall. Ihe parameters P] through P 3 and Qi 
through Q 7 are determined empirically by means of least squares curve fitting 
of equation (18a) with the test data. 


Unsteady Moment Coefficient 


c 


mc/4 



Cue (<* - Aa„) *■ 3 Aa„ + AC 
MS ' 2 om 2 m 


(19a) 


AC 


m 


- n l A + n 2 a w f n 3 (a^ ) f n 4 | a w| + n 5*l * n 6 Aa 2 f n 7 a Dm A 


Dm S m 


(19b) 


In equation (19a), a om is the static moment at zero angle of attack, normally 
zero, and is the static moment coefficient about quarter chord. The 

parameters n-| through are determined by least square curve fitting of 

the test data. For unstalled airfoils, the last three terms are zero. 


It should be noted that ( 1 ) the delay, the time from the onset of static 
stall to actual stall, is taken into account by calculating the actual dynamic 
stall angle from the curve fitted experimental data, ( 2 ) the difference in 
pitching and plunging motion is not taken into account explicitly. Ihe total 
angle of attack and its derivative are used in the formulation, and (3) vortex 
formation and its effects are accounted for. 


MODEL C ONERA MODEL 

The ONERA model (ref. 39) describes the lift and moment coefficients in 
terms of ordinary differential equations. The equations are of first order in 
the linear flow regime (attached flow), and they are of third order in the 
nonlinear regime (separated flow), the third order equation simulates the 
pseudo inertial, pseudo damping, and pseudo elastic forces of a fluid oscilla 
tor. The parameters (six for lift, and five for moment) in the differential 
equations are identified, as functions of angle of attack, by parameter iden 
tification of test data. These tests are conducted about each mean angle of 
attack at various reduced frequencies at small amplitude (« 1 °) of oscillation. 
They showed good correlation for angles of attack up to 23°. Reynold's number 
and Mach number effects are included implicitly from a static airfoil data 
base. It should be noted that a differential equation automatically accounts 
for the time-history effects. 

The model consists of three equations that relate the lift coefficient of 
an airfoil to its angle of attack. They are: 

^ ^ r A* A 

C| 1 t \C| 1 - \C| ^ \S 0 r act + SO (20a) 

C ( 2 f aC | 2 * rC l2 = -(r AC j f Ea) (20b) 
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a = 


(20c) 


c 



+ C 


L2 * 


e 


+ 


* 

h 

b 


where C |_i and Cl 2 are the lift coefficients in the linear and nonlinear 
regions of angle of attack (below and above static stall angle respectively), 
a is the total aerodynamic angle of^attack of the airfoil, e is the angle of 
attack due to pitching motion, and h/b is the angle of attack due to 
plunging motion. C| a is the static lift coefficient in the linear region of 
angle of attack. ACl is the difference between the extended linear lift curve 
( C lsl = a 0 8.a) an d the actual static lift curve (Cls). c 8. is the resulting 
total lift coefficient as shown in figure 5. 

Similarly, the moment coefficient can be expressed as a first order dif- 
ferential equation for angles below static stall angle with an additional sec- 
ond order differential equation for angles above static stall angle. However, 
the first order differential equation is not necessary for moment coefficient 
since the Theodorsen function C(k) does not appear in equation for moment 
coefficient in Theodorsen's moment equation for a flat plate, reference 52, 
equation (5-312). Therefore the equations for moment coefficient are given as 


* * ** 


C Ml = C Mst f s0 f oa f se 

(20d) 

¥ 3 ^M2 f r ^M2 = "( r + 

(20e) 

c mc/4 = C M1 * C M2 

(20f) 


where C^] and are the moment coefficients in the linear and nonlinear 

regions of angle of attack. Cmj is the static moment coefficient in the lin- 
ear region. ACj^ is the difference between the extended linear curve and the 
actual static curve as shown in figure 5. c mc //| is the resulting total moment 
coefficient. 

In equations (20a) to (20e), the (*), (**) operators represent deriva 
tives with respect to nondimensional time, t = V t/b. The values of the 
parameters, X, s, a, a, r, E are different for lift and moment coefficients, 
and are obtained from wind tunnel tests. However, their general behavior can 
be identified and is presented below. 

for an airfoil with a fixed Reynold's number, and Mach number, the param- 
eters X, s, a, a, r, E are functions of the blade angle of attack only, and 

must be determined from wind tunnel test data of the corresponding airfoils by 
parameter identification. The parameter, X, is the time-delay parameter asso- 
ciated with the lift deficiency function. It provides for changes in magni- 
tude and phase of the lift. The parameter s is the apparent mass term. The 
parameters in equations (20b) and (20e) are associated with the stall phenome- 
non. In particular, a is a damping parameter; r is the frequency of the 

stall response, and E is a phase shift parameter associated with the stall 

response. The parameters defining the differential equations showed a remark 
ably common behavior with several airfoils, that is they have shown approxi- 
mately the same values for r, a, E, a for each airfoil, and same dependence 
on AC l - The parameter r happens to be the same for lift and for moment. 

This means that the resonance frequency is in fact a property of the flow 
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which manifests Itself In every type of force. The static values (C|_s, Cm$) 
can be approximated by polynomials. The static curves are very sensitive to 
airfoil shape and conditions of flow, but the unsteady characteristics are not. 

It should be noted that the stall delay time, the time from the Instant 
the angle of attack exceeds the static stall angle to the occurrence of actual 
stall, is introduced explicitly while solving the equations. This Is done by 
keeping the right-hand side of equations (20b) and (20e) equal to zero until 
the static stall angle Is reached and the delay time Is met. This delay time 
(In terms of reduced time), based on experimental data, Is taken as 10 for 
lift and 5.8 for moment (see ref. 40). The effects of pitch and plunging 
motion are distinguished by separating the terms that multiply the total angle 
of attack (pitch + plunge) and that multiply the pitch rate terms. This Is In 
contrast to model A and model B, where such distinction was not made. However, 
it should be noted that In the ONERA model, the system dimension Increases with 
number of integration points (3 for lift and 2 for moment), and it can be used 
In a linearized stability analysis. 

In reference 53, the equations (20a) to (20e) are used to obtain the 
dynamic response of a typical section helicopter blade, and In reference 54 to 
an entire helicopter blade. It should be noted that this type of modeling, 
expressing aerodynamic force in the form of differential equations has been 
attempted In reference 55. In particular reference 56 has used It In the 
dynamic stall analysis of helicopter blades. 

A flow chart of the three dynamic stall models Is given in figure 6. 


GOVERNING AE ROE LAST 1C EQUATIONS OF MOTION 

This section presents the governing equations used in the aeroelastic 
study conducted with the three dynamic stall models described earlier. Three 
simple structural models are considered. 

1. Calculation of lift and moment coefficients for airfoils oscillating 
sinusoidally in pitch about quarter chord. This study will show the degree of 
correlation with published results. 

2. A typical section model having pitching and plunging motion. 

3. A plate model having pitching and plunging motion. 

The flat plate case simulates the case of an advanced turboprop while reducing 
the complexity in the structural modeling. When the aeroelastic equations are 
formulated in terms of normal modes, for example as in reference 2, the struc- 
tural modeling is Identical with that of flat plate. For the aerodynamic force 
calculations, the dynamic stall models require as Input the angle of attack 
and Its derivatives at each section along the radius. This can be easily 
Implemented in the flat plate model by Including the rotational effects In the 
calculation of angle of attack and its derivatives. Hence a flat plate case 
is a natural example to start with and to extend to the advanced turboprop. 
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Airfoils Oscillating Sinusoidally in Pitch About Quarter Chord 

For stall models A and C, the angle of attack and its time derivatives 
are given by 


a = a„ + a sin ut 
0 

a = (>> a COS ut 
a = - a Sin ut 


( 21 ) 


where a is the instantaneous angle of attack, <s 0 , mean angle of attack, a 
is amplitude of oscillation, and u is the frequency of oscillation. A time 
delay of 10 for lift and 5.8 for moment is used in model C, based on the study 
of reference 40. 

For stall model B, the angle of attack, a, nondimens ional time derivative 
of angle of attack, A, and the decay parameter, a w , required in the calculation 
of the effective angle of attack are given by 


+ a sin 

(ut) 

(22) 

+ a sin 

(ks) 

(23) 

a COS 

(ks) 

(24) 


where k = u b/V, 

Y-| (k.M) 

Y 2 (k,M) 


a w = ( k , M) f (a - 

s = Vt/b and 

0.165 (1-M 2 ) (0.0455) 

2 2 2 2 

k t (1-Frr (0.0455) 

0. 165 k 2 

2 2 2 2 

k f (1 M ) (0.0455) 


* 0 ) ^ 2 (k,M) 

0,335 (1-H j (0,3) 

f 2 2 2 2 

k r (1-M ) 0.3 

0.33 5 k? 

f ? 2 2 2 

k + (1-M ) 0.3 


(25) 


(26) 


(27) 


For each time step n, the angle of attack, A and a w are calculated. 
Ihe angle of attack is compared with the static stall angle. If this angle is 
above the static stall angle, the lift and moment coefficients are obtained 
from corresponding equations of the three dynamic stall models. Then the 
normal force coefficient, Cfj, is given by c t cos (a). It should be noted 
that in contrast to models A and B, model C requires integration of differen 
tial equations. 


A Typical Section Model 

Ihe typical section is a representative section, usually taken at 
75 percent span. Figure 7, shows a typical section oscillating in pitch and 
plunge. The plunge motion h is positive downwards, and the pitch motion a 
is positive nose up. The equations of motion can be written as 
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mh + Sa + C^h * K^h = Qf, 


(28) 


Sh + I a a + C a a + K a a = Q a (29) 

where m is the mass, S is static balance, is the plunge spring coeffi- 
cient, K a is the pitching spring coefficient, C^, and C a are the mechanical 
damping in plunge and pitch respectively, and Q a are the force and moment 
on the ai rfoi 1 . 


Letting t = cot, and £ = h/b, the equations can be written as 



(30) 


(31) 


With zero mechanical damping, and expressing and Q a in terms of 
normal force, C^, and moment coefficients, c m , respectively, equations (30) 
and 31 can be written as 



* 2 

(u kr 


[k! U- 


-C. 


iryk* 


m 


where k = <ob/V, U* = V/bco a , y = m/irpb 2 . 
The elements of [M] and [K] are 
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= 1 . 0 , 
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= M 
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= x 




(32) 


where x a is the static unbalance, r a is radius of gyration, co^, co a are 
uncoupled bending and torsion frequencies, p is the density of air, b is 
semi chord, m is the mass per unit length, y is the mass ratio, V is the 
resultant velocity, and to is reference frequency. 


In the previous section, the angle of attack and its derivatives were 
given explicitly for sinusoidal oscillations. However, a structure performs a 
more general motion, lherefore, it is required to calculate the angle of 
attack (a), the nond imensional time rate (A) and the decay parameter (a w ), 
for an arbitrary motion. This is given below. 

At any time step n, the velocities tangential and perpendicular to the 
leading edge are given by (see fig. 7) 

Upsh+Xa+V sin (a 0 r a) (33) 
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U-| = V COS (a 0 + a) 


(34) 


from which the instantaneous angle of attack is given by 



(35a) 


With this instantaneous angle of attack, the parameter A and a w are deter- 
mined by using the following equations: 


(A), 


, , . (1.5a - 2a . t 0.5a _) 

( d_a) ' n n-1 n-2 


(ds), 


(As), 


n ' ' n 

using the backward difference scheme, and 

(a w ) n = f ^n 


(35b) 


(35c) 


where 


0.0455 (1 M ) (as) 


X = X , e 
n n- 1 


* °‘ 165 ( *n - Vl> 


-0.3 (1-M ) (As) 


Y - Y , e 
n n-1 


f 0.335 (a p - a n -| ) 


(35d) 


(As) = 2 U ” 

' n n c 

where c is the chord, U n is the velocity normal to the leading edge at 
time step n, and subscripts n, (n 1), etc., indicate the values at the n, 
(n-1) th step in a step by step integration method. 


For these aerodynamic parameters, a, A, a w , the lift and moment coeffi 
cients (about quarter chord) are calculated from the dynamic stall models A 
and B described earlier. Then the aerodynamic normal force and moment coeffi- 
cients about the elastic axis are given by 


C N ■= c* COS a 


c = c 
m mc/4 



0.5) 


(36) 

(37) 


where c mc //| is the moment coefficient calculated about quarter chord 
point. In the case of dynamic stall model C, the stall parameters appearing 
in equations (20a) to (20e) are calculated using the above aerodynamic parame 
ters, and the differential equations (20a) to (20e) are integrated along with 
the structural equations given in equation (32). A state vector form of this 
combined structural/stall model is given in section A.l of appendix A. 
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Plate Model 


lhe governing equations of motion are formulated In terms of normal modes 
of a cantilevered plate. The nonrotating normal modes of the plate are 
obtained from a NASIRAN analysis. Since any general vibratory motion of the 
plate can be expressed In terms of translation of some reference axis and 
rotation about this reference axis, the plunge and pitch motions are expressed 
about a reference axis. If h(y) Is the bending deflection of the reference 
axis at station y, when the blade Is vibrating In the 1th normal (coupled) 
mode, and a(y) is the rotation about the reference axis, the displacement at 
station y can be expressed as a superposition of the contributions of the 
various normal modes as 


h * X>i Q^t) 
a = q i( *•) 


(38) 


where q-j(t) Is the generalized coordinate, which Is a function of time, lhe 
amplitude q^(t) expresses how much of each normal mode Is Introduced Into the 
general vibratory motion. The kinetic energy dT of an element dY is 

dl = \ [m(y)h(y,t) 2 f I (y)d(y,t) 2 ♦ S (y)h(y,t)d(y,t) ] (39) 

C Cl 

Substituting for h and a and integrating over the span the total kinetic 
energy is given by 


where 


1 



.2 




2 Sh.a.) dy 


with the normality condition 



f a a 1 a j f S (^1 a j 


The potential energy is 


a i h j ) ] dy ^ 0 


(40) 


('ll) 


(42) 


U 



CJ . 


(43) 


where w-j is the circular frequency of vibration In the ith normal mode. 
Then the governing equations of motion can be written as 


m i q i 


m.. 




(44) 


where is the generalized force including aerodynamic forces. 
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Aerodynamic Velocity Expressions 

In the present study, simple strip theory is used to calculate the blade 
velocities and forces. Therefore the plate is divided into a number of strips 
(segments), for a flat-plate with no pretwist, the velocities at each segment 
are given by the same expressions as for a typical section model. Therefore 
the angle of attack is given by equation (33), after expressing h and a 
motions in terms of normal modes. Ihen c t and c mc /4 are calculated from 
any of the dynamic stall models described earlier. The lift and moment coef- 
ficients at this segment are given by equation (36), from which the actual 
lift and moment are calculated. Then the generalized force at segment, j, in 
ith mode is given by 


f i j = (hi (y)L ♦ c»i ( y ) M) dl (46) 

where dl is the length of the segment, L and M are the lift and moment at 
jth segment, and hi (y) and ai(y) are values of h and a at the center 
of the jth segment. 

lhe total generalized force on the blade for ith mode is given by 

N 

Q. = T f.. (47) 

1 j“l 1J 

where N is the number of segments, the structure is divided into, lhe 
equations of motion are solved for this time step. 

It should be noted that for stall model C, the lift and moment coeffi 
cients are given in the form of differential equations, and the combined 
structural/stall model has to be solved for each time step. For a plate (con- 
tinuum) model, the number of equations in the combined structural/stall model 
depend not only on the number of normal modes used, but also on the number of 
segments into which the blade is divided. A state vector form of this com 
bined structural/stall model for model C is given in section A. 2 of appendix B. 


RtSULIS AND DISCUSSION 

lhe numerical results of the study are presented in this section, lhe 
three aerodynamic stall models explained earlier are used for comparative 
studies. The results and discussion are performed for three simple structural 
models: (1) an isolated airfoil oscillating sinusoidally in pitch about 

quarter chord; (2) a typical section model performing pitching and plunging 
oscillations; and (3) a flat plate model performing pitching and plunging 
oscillations. A computer program with these dynamic stall models and struc- 
tural models has been developed and will be expanded to consider the aero 
elastic response of advanced turboprops, lhe following results show the 
validity of the computer program and an assessment of the three dynamic stall 
models. The experimental data presented in this report was obtained from 
enlarged published plots. 

In this study all the calculations are made for three airfoils, namely 
NACA 0012, NACA 0012 (MOD) and 0A212, and for Mach numbers of 0.3 and 0.4, 
since the semi empirical data is readily available for these Mach numbers. As 
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pointed out earlier, the dynamic stall effects vary with airfoil geometry and 
Mach number. However, the following comment can be made for stall at higher 
Mach numbers. The static stall angle decreases as the Mach number increases 
and the blade may stall earlier than at low Mach numbers. Experimental data 
of reference 56 showed that for Mach number equal to 0.6 there is shock- 
induced separation and stall. The dynamic data suggested that the formation 
of shock waves somehow inhibit the development of the vortex shedding process 

At this pint, a mention of the effects of other parameters on dynamic 
stall is in order. Experimental results of reference 7 indicated that, (1) 
increasing reduced frequency decreases the intensity of the dynamic stall 
vortex shedding and delays the formation and growth of the leading edge 
vortex, (2) cambering the leading edge tends to delay stall onset, (3) the 
higher the Reynold's number, the later the separation appears and develops, 
(4) sweep reduces dynamic stall effects, (5) the effect of decrease in thick- 
ness is similar to the effect of an increase in Mach number, and (6) cascades 
reduce dynamic stall effect. The dynamic stall effects also depend on mean 
angle and amplitude of oscillation. 

It should be noted that when comparing with published results some quan 
titative differences may exist since the dynamic stall loads are sensitive to 
airfoil static data used. For completeness, the static data (c a and c mc /4 

versus a), and the dynamic stall data for the airfoils mentioned above are 

given in appendices B and C. When comparing results from model C with those 

obtained from model A and model B, it should be noted that model C involves 

integration of lift and moment equation simultaneously with structural 
equations . 


Isolated Airfoil Oscillating in Pitch About Quarter Chord 

Verification of the published resu lts . - Ihe intention here is to verify 
the published results and to check whether the dynamic stall models are cor- 
rectly implemented in the computer program. 

Stall model A: In reference 35, this model was applied to four airfoils, 

namely the V230I0 1.58, NACA 00I2(M00), V13006 .7, and NACA 0006, and the 
theoretical and measured lift and moment coefficients were correlated. In the 
correlation, a variation in the mean angle of attack, amplitude of oscillation, 
Mach numbers, and reduced frequency was considered. The test data included 
forced pitch oscillations for the four airfoils and plunging oscillations for 
the V23010 1.58 airfoil. It was concluded that overall correlation between 
theory and test was good for the normal force coefficient and was acceptable 
for pitching moment coefficient. 

lo facilitate interpretation of analytical model A, the lift and moment 
coefficients are calculated herein and are shown in figures 8 and 9 along with 
the measured data. Ihe NACA 00 1 2 ( MOD ) airfoil in pitching motion at 0.4 Mach 
number is considered. Figure 8 presents the plots for a mean angle of 9.93°, 
amplitude of oscillation of 4.65° at a reduced frequency of 0.064 and figure 9 
presents for a mean angle of 12.25°, amplitude of 4.8°, and a reduced fre- 
quency of 0.126. Ihe figures reproduce the corresponding ones presented in 
reference 35. It can be seen from these figures that (1) the lift coefficient 
is predicted well, (2) the maximum moment coefficient has not been obtained 
for the case of k - 0.126, and (3) moment loops are not reproduced correctly. 
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A close examination of the results in reference 35 for other airfoils, and for 
other flow conditions revealed that the prediction is good for thin airfoils 
at low reduced frequency and the maximum moment coefficient is not predicted. 

Stall model B: Reference 38 presents synthesized lift and moment loops 

for 13 data sets comprising di f f erent ai rf oi Is (SC1096, NACA 0012, V0012, VR7, 
NLR-1, V2301 1.58, Yawed 0012), Mach numbers (0.18 to 0.6), reduced f requen 
cies (0.0 to 0.38) and Reynold's numbers 2.5 to 6.2 million). In the present 
report NACA 0012 is chosen to verify the program. Figures 10 and 11 present 
lift and moment loops for two cases. The cases represent different mean 
angles, amplitude, and reduced frequency at 0.3 Mach number and Reynold's num 
ber of 3.8 million. The steady state-data for this airfoil is obtained from 
an analytical curve fit program available in G400PR0P computer code (ref. 57). 
Ihe static stall angle for Mach number 0.3 is 12°. The lift loops from this 
model show good correlation with experimental data over the upper portion 
(upstroke) but the present program under predicts the lift coefficient in the 
returning stroke. But the published results in reference 38 showed good com 
parison even in the down stroke. Ihe reason for this discrepancy between the 
present theoretical results and those of reference 38 is not known at this 
time of writing. However, it should be noted that in a response analysis, one 
is interested in calculating the maximum lift obtained in dynamic stall. So 
the predicted loops are assumed satisfactory for this purpose. Ihe moment 
loops show good correlation with the theoretical and experimental loops of 
reference 38. In contrast to lift loops during the down stroke, the moment 
loops behave well throughout the cycle. The good correlation for lift and 
moment loops may not be surprising since the model is actually experimental 
data, reconstructed by analytical curve fitted expressions, which naturally 
takes into account all the events that occur in a dynamic stall process. 

Stall model C: The ONFRA model was first fitted for the 0A2I2 airfoil 
(ref. 39). This airfoil has a lift stall angle of 10° and a moment stall 
angle of 6°. Ihe results are presented for a Mach number of 0.3 and reduced 
frequency of 0.05, with amplitude of oscillation of 6° for mean angles of 
attack of 12° and M°. Ihe governing equations are solved using a fourth 
order Rurtge Kulta method. Figures 12 and 13 with designation "a" show the 
lift loops obtained from this model with and without stall delay for mean 
angles of 12° arid If. Ihe model with stall delay predicted the maximum lift 
coefficient well. The lift loops calculated herein are in good agreement with 
corresponding ones in reference 39. The moment loops for the same cases are 
presented in figures 12 and 13 with designation "b". Again, the results show 
good correlation with those published in reference 39. Introducing stall 
delay did not improve the correlation since the range of angle of attack 
studied is well above the moment stall angle of 6°. It is interesting to note 
that the moment curves (even the experimental ones) did not show any loops for 
this airfoil. 

Comparison of the three dynamic sta ll models f or a single airfoil. - In 
this section, the three dynamic stall models are applied for a particular air 
foil, lri this case, NACA 0012 airfoil, at Mach number = 0.3, and Reynold's 
number of 3.8 million is considered. 

For NACA 0012 airfoil, the lift loops from model B are presented earlier, 
arid they are calculated here for models A and C. In generating these plots, 
the static stall angle for lift is assumed as 12° for models A and B, and 14° 
for model C. To compare all these results, all three loops are shown in 


figures 14 to 17, for four cases by varying the mean angle of attack, amplitude 
of oscillation, and reduced frequency. These cases are listed below. 


Case 

1 : 

a 0 = 10° 

i = 10°, 

k = 0.04813 

Case 

2: 

a 0 = 12°, 

a = 10°, 

k = 0.09756 

Case 

3: 

a o = 12°, 

a = 8°, 

k = 0.12528 

Case 

4: 

a 0 = 1 5° , 

i = 5°, 

k = 0.15106 


The measured lift loops are also Included In figures 14 to 17. 
obtained from dynamic stall models A, B, and C are designated by a, 
respectively. To facilitate comparison, the experimental plots are 
In the figures. 


The plots 
b, and c, 
also shown 


Figure 14 shows the lift loops obtained for case 1. All three stall 

models predict the loops reasonably well for this reduced frequency of 0.048. 
Even the simplest model, model A, gives a very good correlation between theory 
and experiment. Model B, as mentioned In the previous section underpredicts 
the lift coefficient during the down stroke. The ONERA model, model C, 
predicts a high lift coefficient compared to models A, B, and experiment. 

This may be due to the assumption of high static stall angle for the airfoil, 
and may be due to the numerical procedure required In Integrating the 
equations . 


Figure 15 presents the plots for case 2. The reduced frequency for this 
case Is approximately twice that for the previous case. The amplitude of 
oscillation Is 10°, as In the previous case. It can be seen from plot 'a' 
that model A did not qualitatively reproduce the experimental curve, even 
though the maximum lift coefficient Is reasonably well predicted. Since the 
area of the loop represents the energy transfer, or damping In the system, 
this model underpredicts the damping. Models B and C predict the loop 
reasonably well. 


Figure 16 shows the lift loops obtained for case 3. The reduced frequency 
and the amplitude of oscillation are 0.126 and 8° respectively. Models A and 
C predict a larger maximum lift coefficient and model B predicts a lower maxi 
mum lift coefficient compared to experiment. 


Figure 17 shows the plots for case 4. The reduced frequency Is 0.161, 
and the amplitude of oscillation Is 6°. It Is surprising to note that model A 
predicts the lift loop reasonably well compared to experiment. Model C pre- 
dicts a larger lift loop. 

An overall comparison of the figures 14 to 17 Indicate that all three 
dynamic stall models predicted the lift loops reasonably well for the cases 
considered. However, the following note Is In order regarding the application 
of model C. Ihe parameters in model C have been evaluated from tests con- 
ducted at small amplitude of oscillation (»1°). Here the model has been 
used for predicting lift at high amplitude of oscillation ( >5 0 ) . This model 
is expected to give still better correlation with experimental data when used 
In a small amplitude of oscillation environment, such as that exists for pro 
pellers and compressors. Models A and B do not have such restrictions on 
formulation and use. 
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Moment Coefficient 


Figures 18 to 21 show the moment loops obtained from models A and 8. 
Specifically, in these figures, plots (a) and (b) show the loops obtained from 
models A and B, respectively. Since the moment data is not available for 
NACA 0012 from the ONtRA fit, loops from this model are not shown. Ihe plots 
are presented for mean angles of attack of 10°, 12°, and 15°, amplitude of 
oscillation of 10° and 5°, and reduced frequencies of 0.09766, 0.09633, and 
0.161. An overall observation of the figures 18 to 21 shows that the moment 
loops are not predicted well by model A. Ihe experimental data shows three 
loops. Model A predicts only two loops. Since the area of the curve is pro 
portional to the damping in the system (or energy dissipated/absorbed), stall 
response curves from model A may not be accurate. For model B, there is a 
very good correlation between calculated and measured results, lhis should be 
expected with this model, because these loops are curve fitted to experimental 
data . 


Aoroelastic Response Studies with Dynamic Stall Models 

Ihe following sections present the results of the response study conducted 
with two structural models, a typical section model and a plate model. Both 
models perform plunging and pitching oscillations. For a given airfoil, the 
comparison of the limit cycle behavior with both pitching and plunging oscil 
lations from the three dynamic stall models would have been ideal. However, 
the moment data is not available for NACA 0012 from model C, for a comparative 
study. Iherefore response from plunging only is considered for comparison 
purpose of the three dynamic stall models. However, for models A and B, both 
plunging and pitching are considered and compared. Ihe governing equations 
are solved by the Wilson 0 method (ref. 69). This method assumes a linear 
variation of acceleration between two time steps. Ihe method is an implicit 
integration method and is unconditionally stable. 

Plunging motion only. - Ihe NACA 0012 airfoil at a Mach number of 0.3 is 
considered for the numerical study. Ihe response curves are obtained for zero 
initial conditions with a time step equal to or less than the l/10th of the 
minimum period of oscillation. Ihe variation of the angle of attack with time 
is given by a = a 0 t- h/V, where h is the plunging displacement. 


Typical Section Model 

Ihe governing differential equation for this case is equation (30). Ihe 
values of the parameters in the equation are, mass ratio, p - 76, w = 40.0 
rad/sec, and - 66.9 rad/sec. Figure 22 shows the variation of plunging 
displacement with time at an initial angle of attack of 4.6° obtained from the 
three models. Ihe response from the three models shows' a converging trend 
towards a steady value. Figure 23 shows the variation of plunging displacement 
with time when the initial angle of attack is 16°. Ihe models predict differ- 
ent behavior. Model A shows a small amplitude limit cycle behavior, model B 
shows a converged solution, lhis is due to the fact that model B includes 
unsteady aerodynamic effects, and hence predicts a higher stall flutter speed. 
Model C shows a high amplitude limit cycle behavior since the model in general 
overpredicts the lift for this airfoil as shown in the previous section. 
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Plate Model 


tquatlon (44) is the governing equation of motion for this case with 
i = 1. The values of m-| , and w| are 0.972x10'^ lb sec^/in. and 
999.0 rad/sec. Figures 24 and 25 present the plunging response obtained from 
the plate model. Figure 24 shows the response obtained from the three dynamic 
stall models for an initial angle of attack of 4.5°. All three dynamic stall 
models predict a converging solution. Figure 25 shows the response obtained 
for an initial angle of attack of 15°. The behavior is as that predicted 
above for the typical section model. 


P ltc hi ng and P lungi ng Oscillation s 

Typical Section Model 

The governing equations for this combined bending and torsion are given 
by equation (32). The numerical study is conducted for the following 
structural parameters. 

p = 76.0, x a = 0.25, a h = 0.15, r a = 0.6229, 

up = 55.9 rad/sec, u ft = 64.1 rad/sec, Ch - 0-0. 

C a = 0.0, b = 0.127 m (5") 

For these values, the model has both inertial coupling and aerodynamic cou 
pling between the plunging and pitching motions. Before starting the response 
analysis, it was necessary to find the aeroelastic parameters corresponding to 
a neutrally stable condition. Hence, a classical flutter analysis was per- 
formed in the frequency domain by using two-dimensional unsteady aerodynamic 
theory. This predicted the flutter speed, Vp , as 27.46 m/sec (90.1 ft/sec), 
and reduced frequency, kp , as 0.27624. lhe corresponding velocity parameter. 

Up , is given by 

* V F 

U, -- r - 3.3/35 
F bu 

a 

and this flutter speed corresponds to a Mach number of 0.0811. 

In the present time domain analysis, the response analysis was carried out 
by varying U*. keeping the values of other parameters constant. The static 
data corresponds to a NACA 0012 airfoil, at a Mach number of 0.3 and Reynold's 
number of 3.8 million. The static stall angle for this airfoil is 12°. lhe 
initial conditions are taken as = 0., £ = 0., a 1 = 0., and a = 0.01 rad. 

A time step equal to or less than l/10th of the minimum period of oscillation 
is used in the calculations. 

Figure ?6a shows the variation of pitching displacement with time when 
the airfoil is at an initial angle of 4.5°, obtained from model A. It should 
be noted that for this angle of attack, which is below the static stall angle, 
the aerodynamic theory in model A corresponds to quasi steady theory, lhe 
plot shows three graphs that correspond to U* = 2.50, 2.675 and 2.75. The 
graph corresponding to U* ^ 2.50, shows converging oscillations, whereas 
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the one for U* = 2.675, shows sinusoidal oscillations. This Is assumed to 
be Indicating neutral stability, since the angle of attack variation for this 
case Is well below the static stall angle. The graph for U* = 2.75, shows 
diverging oscillations. The value of U* = 2.675 Is taken as the one cor- 
responding to the flutter speed. This Is 21 percent less than the value pre- 
dicted from unsteady theory, showing that quasi-steady theory underpredicts 
the flutter speed. Figure 26b shows the variation of plunging displacement 
with time. It exhibits the same behavior as the pitching displacement. 

Figure 27 shows the response curves obtained from model B for the same 
structural parameters and for the same Initial angle of attack. It should be 
noted that the unsteady parameters (below stall) available for NACA 0012 at a 
Mach number of 0.3 are used In the response calculation. Figure 27a shows two 
graphs that correspond to U* = 3.17, and 3.18. The graph corresponding to 
U* = 3.17, shows converging oscillations, whereas the one for U* = 3.18, 
shows diverging oscillations Indicating Instability of the motion. Flence the 
value of U* = 3.175 corresponds to flutter speed. This Is 6 percent less 
than the value predicted from unsteady theory In the frequency domain. 

Figure 27b shows the variation of plunging displacement with time for the same 
parameters as above. It exhibits the same behavior as the pitching 
displacement. 

Figures 28 and 29 show the response curves obtained from models A and B 
for an Initial angle of attack of 15°. This Is above the static stall angle 
of 12° for NACA 0012. The other parameters used In stalled airloads calcula- 
tion are Re ^ 3.8 million and Mach number = 0.3. 

Figure 28 shows the variation of pitching displacement with time when the 
airfoil Is at an Initial angle of 15°. Figure 28a Is the response curve 
obtained from model A. The plot shows two graphs that correspond to U* 
equal to 0.1 and 0.2. Both the graphs show that the transients (due to non 
zero Initial acceleration and nonzero Initial displacement) die out In a 
certain amount of time and then the oscillation starts growing Indicating 
unstable oscillations. It Is seen from figure 28a that the time for the 
transients to die out is inversely proportional to the velocity parameter, 

U*. Response curves generated with zero initial conditions (zero velocity 
and zero displacement), not shown here, Indicated that the oscillations start 
to diverge immediately, lhis indicates that flutter begins essentially at 
zero speed compared to U* = 2.75 at an initial angle of 4.5°. Flowever, 
this physically unreasonable result is not surprising, since this model did 
not predict the moment loops correctly for this airfoil which resulted in 
negative damping. 

Figure 28b is the response curve obtained from model B, for the same 
structural and aerodynamic parameters. The plot shows two graphs, one for 
U* = 0.45 and the other for U* = 0.5. The graph for U* = 0.45 shows decay- 
ing oscillations Indicating stable oscillations. Ihe graph for U* = 0.5 shows 
that the airfoil Is set into a limit cycle behavior. The decay of the transi- 
ent and the motion being set Into a limit cycle is more clearly seen here, 
compared to that shown in figure 28a obtained from model A. The critical 
stall flutter speed corresponds to U* = 0.475, the average of the above two 
values. This Is 85 percent reduction in flutter speed compared to 3.175 at an 
initial angle of 4.5°. 
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Figure 29 shows the variation of plunging displacement with time, and it 
shows the same variation with time as pitching motion for both the models. 


Plate Model 

Ihe response of the plate is analyzed using normal modes. These normal 
modes are obtained from a NASI RAN analysis. The plate is 3 in. in length, 

1 in. in width (chord) and 0.05 in. thick. The material density is 
2.6xl0' 4 lb sec^/in 4 . Young's modulus 7.89x10^ lb/in? and Poisson's ratio 0.3. 
Ihe analysis used QUAD4 elements. In the present study, two normal modes are 
used in the analysis. For this problem, the first and second normal modes are 
essentially the first bending and the first torsion modes. From the NASIRAN 
analysis, the first two generalized masses are 0.9720xl0' 5 lb sec 2 /in. and 
0.161x10 ^ lb sec 2 /in., and the corresponding frequencies are 999.2 rad/sec 
and 5860.55 rad/sec. The NACA 0012 airfoil at Mach number 0.3 is considered 
for the aeroelastic study. Ihe governing equations of motion are given by 
equation 44. Ihe response is calculated with the following initial conditions: 

I I 

q.j = q.| = q^ = 0.0 and q^ = 0.01, where q^ and q^, represent the first two 

I I 

normal coordinates and q^ and q^ their derivative with nondimensional time. 
The variation of the normal coordinates with time is shown in figures 30 to 33. 

Figure 30 shows the variation of the first normal coordinate (the first 
bending mode) with time when the flow is approaching the plate at an initial 
angle of attack of 4.5°. Since the initial angle of attack is below static 
stall angle (12° for this airfoil), the response is a quasi-steady response 
for model A. Both models predict the same qualitative behavior, even though 
model B predicts a slightly lower amplitude due to unsteady effects included 
in the model. Ihe response is seen to decay indicating an approach to steady 

displacement. Figure 31 shows the variation of the second normal coordinate 

(first torsion) with time at the same initial angle of attack. Ihe response 
again seen converging to a steady value. 

Figures 32 and 33 show the variation of first and second normal coordi 
nates, respectively, with time when the initial angle of attack is 15°, which 
is above the static stall angle. Both models predict a diverging type of 
oscillation. This means that the plate is unstable for this initial angle of 
attack. Further investigation is required to predict the stall flutter bound 
ary from both models with other initial angles of attack varying from 4.5° to 
15°. Flowcver, it can be noted from previous sections that model B predicts a 

higher initial angle of attack than model A lor the plate to become unstable. 


DISCUSSION OF ADVANCLU 1URB0PR0P APPLICA I ION 

Ihe theory and the computer program presented in this report can be 
directly applied to predict the response of an advanced turboprop with dynamic 
stall. However, proper blade data for the advanced turboprop airfoil is 
needed. Figure 34 shows the airfoil section distribution for SR2 advanced 
turboprop. Ihe thickness ratio ( t/b) , twist (All), design lift coefficient 
( C | □ ) , and planform (b/D) distribution were established to provide for high 
efficiency. Ihe airfoil sections used are NACA 16 series from tip to the 
45 percent radius portion, and 65 series with circular arc (CA) camber lines 


29 


from 37 percent radius to the root with a transition fairing in between. 

These airfoils have been chosen for their high critical Mach number and wide, 
low drag buckets. This necessitates generation of static and dynamic data for 
a number of airfoils along the radius. 


lhe operating conditions of the advanced turboprop indicate that the 
advanced turboprop blades are subject to high frequency, low amplitude oscil- 
lations. This is a favorable condition for light dynamic stall, since the 
penetration of the airfoil into stall is small, and the duration in the stall 
condition is small, even when it is operating at static stall angles. This 
does not allow the boundary layer to reach the leading edge or the formation 
of a vortex like disturbance. This implies that models A and C may predict 
the stall behavior of advanced turboprop reasonably well. Also as mentioned 
in earlier sections, model C has been formulated under the assumption of small 
amplitude, from the theoretical analysis point of view, boundary layer meth 
ods, rather than full Navier- Stokes solutions may be sufficient for advanced 
turboprop application. However, it should be cautioned that light dynamic 
stall is very sensitive, to airfoil shape, particularly, and to the leading 
edge radius. Depending on this radius, the flow may separate at the leading 
edge and may result in deep stall. It is worthwhile here to mention the con- 
clusions drawn from an experimental study of stall flutter characteristics of 
NACA 16 series in reference 60. The study indicated that: (1) an increase in 

camber, thickness ratio or design lift coefficient increases the static stall 
angle, (2) an increase in leading edge radius or in the droop has a slightly 
detrimental effect on stall, and (3) an addition of trailing edge camber 
increases normal force coefficient before stall. 


Ari extensive literature search has found that no static data at and 
beyond stall and no dynamic data exists for 16 series airfoil. Hence the 
following procedure is recommended to obtain static and the dynamic data to 
use for advanced turboprop analysis. Generate the steady data and a single 
dynamic loop data for a representative NACA 16 series airfoil using Navier 
Stokes solution methods (ref. 17). lhen develop the ONtRA dynamic stall model 
from the data obtained from the Navier Stokes code. Use this ONtRA Navier 
Stokes model to calculate the stall induced loads at all 
radius when the angle of attack exceeds the static stall 
the response ol the advanced turboprop obtained from the 
dynamic stall model, and model A. 


sections along the 
angle. Then compare 
ONtRA Navier Stokes 


SUMMARY AND CONCLUSIONS 

Three dynamic stall models are compared for performance and easy imple- 
mentation. The response of a typical section model and a plate model are 
studied for these stall models. The procedure and the computer program can be 
directly used to analyze the response of advanced turboprops with dynamic 
stall. 

In summary: 

1. All three dynamic stall models are implemented in a computer program 
and the program is checked out. 
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2. All three dynamic stall models are applied to a single airfoil to cal- 
culate lift and where possible moment loops, and the strong and weak points of 
these models in predicting the stall behavior are Identified and discussed. 

3. Stall flutter response is calculated by using all three dynamic stall 
models and two structural models. 

4. A formulation for the stall flutter analysis code for advanced turbo 
prop is suggested. 

The following conclusions are drawn from this study: 

1. The operating environment of an advanced turboprop favors the condi 
tion of light stall, where the vortex induced loads are not very severe. 
Therefore, simple dynamic stall models like model A and model C can be used 
for advanced turboprop application. 

2. The ONERA model, model C, involves fewer parameters in modeling the 
dynamic stall model than the synthesis procedure, model B. Model A involves 
only one empirical parameter but must be tested for advanced turboprop air 
foils and high frequency low amplitude oscillations that occur in advanced 
turboprops . 

3. The differential form of the ONERA model suggests that it could be 
used in a linearized stability analysis. 


APPENDIX A 


AEROELAST 1C EQUATIONS WITH ONERA DYNAMIC STALL MODEL (MODEL C) 

lhis appendix presents the state vector form of the combined structural- 
stall model for a typical section model and a plate model with ONERA dynamic 
stall model (model C). 


A.l Typical Section Model 

Ihe governing equations of the typical section are given by equation (32), 
and can be written as (with multiplication factors absorbed into the 
matrices) + 


u 

II 

E 

■ 

E 


■ COS a 

0.0 

CM] 

- [C] 


«■ CK] 





( a p 


a P 


a P 


(a h ► 0.5)cos a 

2.0 



(A. 1.1) 


where a is the angle of attack, and operator ()' indicates derivative 
with respect to t (= ut). 

From ONERA formulation, c^ and c mt //| are given by equations (20c) 
and (20f) as 



(A. 1.2) 


Substituting equation (A. 1.2) in equation (A.T.l), we obtain 


f Note 1: a p is the pitching degree of freedom, subscript p is used to 

distinguish it from a, the instantaneous angle of attack. 

Note 2: lhe transformation from the operator ()* to ()' is ()* - k( ) ' , 

where k = « = reference frequency, b is semi- chord and V is reference 

veloci ty . 

Note 3: lhe multiplication factors include 0.5 pV 2 c and 0.5 pV 2 c 2 to get 

actual lift and moment. 
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B b 3 






(a. f 0.5)cos a 
h 
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0.0 

0.0 

0.0 

1 .0 

0.0_ 



(A. I .3) 


(A.l .4) 


where C, , , C, 2 . c[ 2 , C M| . C^, are given by equations (20a) to (201) as 

' X i X | C | s. 1 ' , '' 

C U 1 ic C Li ' k - Vl° ’ V ‘ S L k0 


•• a L 1 r L 

C | 2 * k C L2 k 2 C l. 2 


k \ AC I. f k L *' 


2 


C M1 * C MSL ^ S M k ° f °M“ k * S M k ~° 
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C uo ♦ r - C,.., «- ~ C..„ = - k‘r.. AC., f , n ~ 


'M2 k M2 k 2 M2 


M M k 


( A . 1 .5) 


where subscripts L and M are used to indicate lift and moment coefficients. 

Combining equations (A. 1.4) and (A. 1.5) and writing it in state vector 
form, we obtain 
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(A. 1.6) 


Del i n i ng 


V -- jf, <* p . K . oy C L1 , C |? , C |2 , C H2 , C M 2 j, 

equation (A. 1.6) can be written as 

(V ) 1 | A ] { V } . IG} 


(A. I ./) 
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A. 2 Plate Model (Continuum Model) 


The governing equations of motion, formulated in terms of normal coordi 
nates, are given by equation 44 as 


2 

ro.q. rrho). q . - Q. 


i = 1,2, ... NM 


(A. 2.1) 


where m^ , , q^ , are respectively the generalized mass, circular fre- 

quency, normal coordinate, generalized force in 1*h normal mode, and NM is 
the number of normal modes used in the analysis. 

The generalized force, , is given by the following equation 

Q i = /* l L(y,t) h . ( y)dy * /* l M(y,t) a. (y)dy (A. 2. 2) 

Jo Jo 

where L(y,t) and M(y,t) are the normal force and pitching moment at y at 
time t, 2. is the length of the plate, h^ and a} are the normal modes in 
plunging and pitching in ith mode. 

Let the plate be divided into N number of segments. Then can be 

evaluated as 


N 


Q, = E [h ( (,.)d,.l..(» r t) . <. 1 (y J )dy J « J (y r t)] 

where j indicates the jth segment, dy j is the length of the 
segment, h -^ ( y j ) and a-j(yi) are the normal modal values at station 
y j . Equation (A. 2.1) for NM normal modes can be written as 


(A. 2. 3) 


[ M ] { q } ♦ [K J {q} = {Q } 


where 


I - 2 2 

[K ] = m-jO^ , m^w^ , 


2 

m.u. , 

i -] 


[M] = Hri-j , , . . . m. , ... » 

{q } = (q r q 2 , ... q y ... q NM ) and {Q} = (Q 1 , Q 2 , ... Q i 


• m MN u NMj 

•• w 
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The vector {Q} Is given by 


Number of segments 


k ' 

- 

O) 

■O 

o 

h 1 (y 1 )dy 1 

t» 1 (y 1 )dy 1 

h 1 (y 2 )dy 2 
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E 
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M j(Vj • t) 1 


t N (y N ,t) 

M N (y N ,t) 

I 

2N * 1 


or 


IQ} - [ 6 6 ]{F} 


( A . 2 . 4 . 2 ) 


I he lift and moment at each segment are related to their coefficients as 
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{F } - { A c p } f [ A 5 p ] {X} 


The coefficients 
given by 



/ 


c 8.1* c ml » ••• c s,j • c mj of each segment are 





L 


1 








(A. 2. 5) 


(A. 2. 6) 


3 / 


or 


{X} = [B 5P J {Z} 

where Z = {Z^ ... 1 . , ... Z N } and {Z.} = |c u , c |2 , c, 2 , c M2 , c M2 J ^ 

Here A,. and B, are the A. and B, matrices of j ^ segment defined 

earlier in Section A.l. Ihe vector Zj is given by equation A. 1.5 for each 

segment. Substituting from equations A. l.b, A. 2. 6, A. 2. 5, A. 2. 4 into 
equation A. 2.1, and writing in state vector form, the final equations are 
given by A. 2. 7, given below. 
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(A. 2. 7) 


where k] , k 2 , ... k j , 
each segment. The size of 


kfj are the values of the reduced frequency for 
the matrix Is ( 2 x NM + 5N). 
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APPENDIX B 


DYNAMIC STALL DATA FOR MODEL B AND MODEL C 


This appendix presents the dynamic stall data for model B and model C 
that Is used In generating the hysteresis loops for the NACA 0012 airfoil. 
The Mach number Is 0.3 and the Reynold's number Is 3.8xl0&. 

Model B: The following data corresponds to data set number 2 In 

reference 38 generated through synthesl zatlon method. The values of the 
parameters appearing In equations 16 to 19 are given below: 


a ss 


12 ° 


Pi /rad = 17.988, 


P 2 /rad - 2.410, P 3 = -0.1661 


Ql /rad = 1.8347, Q 2 /rad = 0.2123, Q 3 = -0.0307, 
Q 5 = 0.7/54, Qfc/rad = 2.842, Q 7 /rad/rad = 7.55/ 


n-j /rad = -1 .782, 


n 2 /rad = -0.191, 


= -0.0035, 


n, = 0.3462, 

D 


n./rad = -1 .405, 
b 


n^/rad/rad = -6.354 


Q/l » -0.0966, 
n^/rad = 0.0082, 


c = 0.125, C Am = 3.64, C wm = -0.05, C At = 0.084, C at = 0.0073, 
C AR = 1.790, C wR - -0.743. 


where ( ) = ( ) a ss . 

Model C: Reference 58 generated the lift hysteresis loop for NACA 0012 

airfoil using a single experimental loop, lhe values of the parameters 
appearing in equation 20 are given as 

X = 0.2, s = 0.09, o = 0.080 0.13AC| 

C| 4 = 0.01 r 0.114 a 

C LS = C L«. » a 5 10° 

C|S = -0.01 «- 0.114 a 0.013 (a - 10 . ) ? - 10° < a < 14° 

C L3 = 1-378 f 0.478 (e" 1 - 5 ^ 14 ) - 1) a > 14° 

Vr = 0.20 «• O.IOAC^ 
a , 0.25 f O.IOAC^ 
t - -0.074CJ, »C L - C u - C |S 

x d (stall delay) - 10, (AC l - 0.208 + 0.0104 (a-14)) 
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APPENDIX C 


AIRE OIL STATIC DA I A 


In this report three airfoils, namely 0A212, NACA 0012, NACA 001 2( MOD ) , 
are used for the aeroelastic response study with dynamic stall. For complete- 
ness, the static data of these airfoils are presented below. The data for 
0A2I2 is presented in polynomial form, and for NACA 0012 and NACA 001 2 ( MOD ) , 
the data is given in tabular form. 


OA212 airfoil: Mach number = 0.3 


Lift coefficient 


ss 


'LH 


'LS 


'LS 


10 ° 

7 - ] 180 * 

7,1 180 a ( a ss < a < a ss } 

7 

£ a. ( | a | 10°) 1 (a < a < 26) 

i-1 

7 

£ a. ( | a | lO 0 ) 1 (-26 < a < a ) 

i =1 


a 0 = 1.24, a! , 0.124, a 2 = 0.063069/, a 3 = 0.01395201, 
a/, -0.0017390851, a b = 0.0001245913, a b = - 4 .6842925x10- 6 , 
a 7 - 7.0879/3x10 8 


For angles greater than 26° (or < 26°), this polynomial diverges. 

Therefore, beyond, |a| > 26°, a constant value of 1.2b is used. 

AC b - 0 . for ( < 01 ^ ^ss) 

= C L9. C LS ( l a I > a ss) 


Moment coefficient 


a SS = 8 ° 

Cmj - L f Pa 

C MS = C MSL * aC M 

where £ = 0.006, p - 0.0014 . I he difference in the linear extension and 

• true static moment, ACf/|, is given below: 

Let O] -- 6 °, 0 2 - 13°, 0 3 - 14.6°, O 4 = 20°, Max = 0.0155 and Min = 0.056. 

Ihen for angle of attack, a, 

between O] and 0 -| ACju) = 0 . 

(a O,)^ 

between 0 | and 0 2 ACm = Max ~ 

( 0 2 - 0 ] ) 
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between o? and 



between o„ and o„ 
3 4 


AC m = -(q - P ) ( a e 3 ) 


> 0 


AC 


Min •• 


0 . 01 ( a - 0 ) 



where 


2Max 


2 


© 


q 


Min t- p(o 4 © 3 ) 
°4 °3 


NACA 00 1 2 ( MOD ) airfoil: The data is presented for two Mach numbers, 

lhe c a and c mc //| with a are given in a, (or c mc /^) pairs. a ss 
is static stall angle, a 0 ] is lift curve slope, and a 0fn is moment curve 
slope, a is in degrees, c& and c mc /4 are per degree. 

lhe following static data is for a Mach number of 0.4. 
a ss - 9°, a 0 n ^ 6.6463/rad, a om - 0. 


(18.86, 1 .042) 

Data c mc/4 /(0.0, 0.), (1.66, 0.), (3.64, -.001), (4.61, .006), (6.6, .008), 

(7.44, .013), (9.34, .063), (11.26, .068), (14.96, -.123), (18.12, -.127) 

( 18.86, .134) , (40. , - .134) 

NACA 0012 airfoil: lhe data is presented for two Mach numbers. The Cj_ 
and c mc /4 with a are given in a, c a (or c mc / 4 ) pairs. a ss is static 
stall angle, a 0 i is lift curve slope, and a om is moment curve slope, a is in 
degrees, c^ and c mc /4 are per degree. 

lhe following static data is for a Mach number of 0.3. 

a ss 3 12°, = 6.30264/rad, aQ m - 0. 



Data c fIIC //) / ( 0 . , 0.), (1.74, - .018), (3.66, .013), (6.66, -.016), 

(7.66, -.001), (9.48, 0.002), (11.37, -.096), (13.36, -.126), (15.16, 

.127) , ( 19.08, .1 52) , (23.94, .182) 


lhe following static data is for a Mach number of 0.6. 


= 5.6°, a o!l = 8.48/rad, a om = 0. 
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Data c 9 / ( 0 . , 0.), (.5, .05695), (1., .11 39), (2.0, .227), (4., .452), (5., 
.563), ( 6 ., .6/2), ( 8 ., .883), (9., .983), (10., 1.0/9), ( 11 ., 1.168), 
(11.5, 1.209), (12., 1.246), (12.5, 1.2/5/), (13., 1.293), (13.5, 1.28), 
(14., 1.229), (14.5, 1.1376), (15., 1.066), (15.5, 1.0343), (16., 1.0097), 
(17., .9255), (19., .9369), (21., .948) 

Data c mc/ 4 /( 0 . , 0.), (1., .00002), (2.0, .00004), (4., .000051), (5., .000032), 
( 6 ., .0), (6.5, -.000014), (/., .000031), ( 8 ., .00005/), (9., .00004/), 

(9.5, .000015), (10, .00004), (10.5, .00012), ( 11 ., .000284), ( 12 ., 

.000164), (12.5, .0005565), (13., -.002985), (13.5, -.00974), (14., 

.02495), (15., .0665), (16., .0/819), (18., -.08931), (20.0, .09839), 

(22. , • .10855) 

Ihe following static data is for a Mach number 0.4. 

a ss - 9°, d o!l 8.944/rad, a om - 0.0015 

Data c a /0., 0.), (.5, .05928), (I., .1185), (2.0, .23652), (4., .46905), 

(5., .5823/), ( 6 ., .692/), ( 8 ., .89991), (9., .99295), (10., 1.0/3/), 

(II., 1.1321), (11.5, I . 1 4 b 4 ) , ( 1 2., 1.1459), ( 12.5, 1.1296), (13., 

1.106), (13.5, 1.09), (14., 1.0865), (14.5, 1.0851), (15., 1.0736). 

(15.5, 1.0433), (16., 1.0121), (17., 1.0112), (19., 1.0095), ( 21 ., 

1.0078), (25., 1.0043) 

Data c mt / 4 /0., 0.), (I., .00012944), (2.0, .0002558/), (4., .0004/534), 

(5;, .0005414), ( 6 ., .0005415/), (6.5, .00050136), (/., .0004223), 

( 8 ., .00009/8), (9., .0005489/), (9.5, .00102), (10., .0015638), (10.5, 

. 0020b93) , (II., .0022224), (12., .00232), (12.5, .0106/9), (13., 

.024558), (13.5, .04069), (14., .05433), (15., .0/1335), (16., 

.0828), (18., .090655), (20.0, .099564), (22., .10953), (25., .1264/) 
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1 ABLE I. - DYNAMIC STALL PREDICTION METHODS 


4 * 


Method for 
predictions 

Refer- 

ence 

Predictive capability 

To stall onset Stall onset 

Post stall 

Motion 

Other factors 

Unsteady pressure 
distribution 

Unsteady boundary 
layer 

lami nar-turbul ence 
transition 

Flow reversal 

Conpressibi lity 

Leading edge 
criterion 

Trailing edge 
separation 

Shock-wave 

interaction 

Acoustic waves 

Vorticity buildup/ 
shedding 

Vorticity 

transport 

Wake modeling 

Induced effects 
of wake 

Subsequent vortex 
shedding 

Reattachment 

Sinusoidal 

Ramp 

Plunge 

Other 

Sweep/3-d i mens i ona 1 

Reynold’s number 

Blade vortex 
interaction 

Roughness, etc. 

Free stream 
turbulence 

Noise 

Navler Stokes 

15 

a 

a 

b 

a 

b 

a 

a 

b 

b 

a 

a 

a 

a 

a 

a 

a 

b 

b 

b 

b 

C 

b 

b 

b 

b 

16 

a 

a 

a 

a 

a 

a 

a 

b 

b 

a 

a 

a 

a 

a 

a 

b 

a 

b 

b 

b 

b 

b 

b 

b 

b 

17 

a 

a 

b 

a 

a 

a 

a 

C 

b 

a 

a 

a 

a 

a 

a 

a 

a 

a 

b 

c 

d 

a 

b 

b 

b 

Discrete 

vortex 

18 

c 

b 

b 

b 

b 

e 

b 

b 

b 

a 

a 

a 

a 

b 

c 

a 

b 

b 

a 

b 

b 

b 

b 

b 

b 

19 

a 

b 

b 

b 

b 

c 

b 

b 

b 

a 

a 

a 

a 

a 

b 

a 

b 

b 

b 

b 

b 

b 

b 

b 

b 

20 

a 

b 

b 

b 

b 

b 

e 

b 

b 

a 

a 

a 

a 

a 

e 

a 

a 

f 

a 

b 

b 

b 

b 

b 

b 

21 

a 

f 

c 

f 

b 

c 

c 

b 

b 

a 

a 

a 

a 

a 

c 

a 

b 

b 

b 

b 

b 

b 

b 

b 

b 

22 

a 

b 

b 

b 

b 

c 

c 

b 

b 

a 

a 

a 

a 

a 

b 

f 

b 

b 

e 

b 

e 

b 

b 

b 

b 

23 

a 

b 

b 

b 

b 

e 

e 

b 

b 

c 

a 

a 

a 

c 

c 

a 

b 

a 

b 

b 

b 

b 

b 
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b 

e 

e 

b 

b 

a 
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a 

a 

a 

b 

b 

b 

b 

c 

b 

b 

b 

b 

b 

b 

Zonal 

25 
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b 
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a Good level of modeling. 

^Not modeled. 

c Approximate level of modeling. 
^Unknown level of modeling. 
e Very approximate level of modeling, 
f Being developed. 
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FIGURE 1. - WIND TUNNEL MODEL OF AN ADVANCED TURBOPROP. 
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FIGURE 2. - THREE TYPES OF STALL. 
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FIGURE 3. - DYNAMIC STALL EVENTS ON NACA 0012 AIRFOIL. 
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THE VARIATION OF T 2 WITH MACH NUMBER (SCHEMATIC). 
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FIGURE 5. - REPRESENTATION OF C Ljg , C L$ , AC l , C m£ , C ms , AND AC M . 
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FIGURE 6. - FLOW CHART OF THE THREE DYNAMIC STALL MODELS. 
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(A) NORMAL FORCE. 



FIGURE 8. - DYNAMIC COEFFICIENT LOOPS FROM MODEL A FOR NACA 0012 
(MOD), a Q = 9.93°, a = 4.68°, K = 0.064, M = 0.4. 
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PITCHING MOMENT COEFFICIENT, C m NORMAL FORCE COEFFICIEANT 



(B) PITCHING MOMENT. 

FIGURE 9. - DYNAMIC COEFFICIENT LOOPS/ROM MODEL A 
FOR NACA 0012 (MOD), a Q = 12.25°, a = A. 8°, 

K * 0.126, M - 0.4. 
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(B) PITCHING MOMENT. 

FIGURE 10. - DYNAMIC COEFFICIENT LOOPS FROM MODEL B 
FOR NACA 0012, a Q = 12°, a = 10°, K = 09756, 

M = 0.3. 
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(B) PITCHING MOMENT. 

FIGURE 11. - DYNAMIC COEFFICIENT LOOPS FROM MODEL B FOR 
NACA 0012, a Q = 15°, a = 5°, K = 0.151, M = 0.3. 
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(B) PITCHING MOMENT. 

FIGURE 12. - DYNAMIC COEFFICIENT LOOPS FROM MODEL C 
FOR 0A212, a Q = 12°, a = 6°, K = 0.05, M = 0.3. 
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PITCHING MOMENT COEFFICIENT, C m NORMAL 


(A) NORMAL FORCE. 



FIGURE 13. - DYNAMIC COEFFICIENT LOOPS FROM MODEL C FOR 
0A212, a 0 = 14°, a = 6°, K = 0.05, M = 0.3. 
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NORMAL FORCE COEFFICIENT 





(C) MODEL C. 

FIGURE 1A. - DYNAMIC NORMAL FORCE COEFFICIENT LOOPS FOR 
NACA 0012, a Q = 10°. a = 10°, K = 0.0A813, M = 0.3. 
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(A) MODEL A. 




FIGURE 15. - DYNAMIC NORMAL FORCE COEFFICIENT LOOPS FOR 
NACA 0012, a Q = 12°, a = 10°, K = 0.09756, M = 0.3. 
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(C) MODEL C. 

FIGURE 16. - DYNAMIC NORMAL_FORCE COEFFICIENT LOOPS 
FOR NACA 0012, a Q = 12°, a = 8°, K = 0.12528, 

M = 0.3. 
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NORMAL FORCE COEFFICIENT 




FIGURE 17. - DYNAMIC NORMAL FORCE COEFFICIENT LOOPS FOR 
NACA 0012# a n = 15°, a= 5°/ K = 0.15106, M = 0.3. 

u CD-86-23722 


64 


PITCHING MOMENT COEFFICIENT, 
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FIGURE 18. - DYNAMIC PITCHING MOMENT COEFFICIENT LOOPS 
FOR NACA 0012/ a Q = 12°, a = 10°, K = 0.9756. 
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PITCHING MOMENT COEFFICIENT, Cm PITCHING MOMENT COEFFICIENT 



(B) MODEL B. 

FIGURE 19. - DYNAMIC PITCHING MOMENT COEFFICIENT LOOPS 
FOR NACA 0012, a Q = 10°, a = 10°, K = 0.09633. 


CO-86- 23703 


66 


PITCHING MOMENT COEFFICIENT 




(B) MODEL B. 

FIGURE 20. - DYNAMIC PITCHING MOMENT COEFFICIENT 
LOOPS FOR NACA 0012, a Q = 15°, a = 5°, K = 0.151. 
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(A) MODEL A. 



FIGURE 21. - DYNAMIC PITCHING MOMENT COEFFICIENT LOOPS 
FOR NACA 0012, a Q = 12°, a = 8°, K = 0.09675. 
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FIGURE 22. - COMPARISON OF PLUNGING DISPLACEMENT FROM THE 
THREE DYNAMIC STALL MODELS, NACA 0012 AIRFOIL, TYPICAL 
SECTION MODEL, a Q = <1.5°, M = 0.3. 
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NORMAL COORDINATE PLUNGING DISPLACEMENT 



FIGURE 23. - COMPARISON OF PLUNGING DISPLACEMENT FROM THE THREE 
DYNAMIC STALL MODELS, NACA 0012 AIRFOIL, TYPICAL SECTION 
MODEL, a Q = 15°, M = 0.3, M= 76. 
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FIGURE 24. - COMPARISON OF PLUNGING DISPLACEMENT FROM THE 
THREE DYNAMIC STALL MODELS, NACA 0012 AIRFOIL, PLATE 
MODEL, a Q = 4.5°, M = 0.3. 
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FIGURE 25. - COMPARISON OF PLUNGING RESPONSE FROM THE THREE 
DYNAMIC STALL MODELS. NACA 0012 AIRFOIL. PLATE MODEL, 

Gq * 15°. M * 0.3 
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PITCHING DISPLACEMENT, RAD 


U* = 2.50 




(B) PLUNGING DISPLACEMENT. 

FIGURE 26. - RESPONSE OBTAINED FROM MODEL A (TYPICAL SECTION) 
°0 ' 
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PITCHING DISPLACEMENT, rad 


U = 3.18 
U* = 3.17 



(A) PITCHING DISPLACEMENT. 



(B) PLUNGING DISPLACEMENT. 

FIGURE 27. - RESPONSE OBTAINED FROM MODEL B (TYPICAL SECTION) 
a 0 - 1.5°. 
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PITCHING DISPLACEMENT. RAD 




TIME IN ut, RAD 
(B) MODEL B. 


FIGURE 28. - PITCHING DISPLACEMENT RESPONSE (TYPICAL 
SECTION) a 0 = 15°. 
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PLUNGING DISPLACEMENT 



(A) MODEL A. 



TIME IN <jt. RAD 
(B) MODEL B. 


FIGURE 29. - PLUNGING DISPLACEMENT RESPONSE (TYPICAL 
SECTION) a 0 = 15°. 
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FIGURE 32. - FIRST NORMAL COORDINATE RESPONSE. PLATE 
MODEL. a Q - IS 0 , M = 0.3. 
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FIGURE 33. - SECOND NORMAL COORDINATE RESPONSE. PLATE 
MODEL, a„ = 15°, M = 0.3. 
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FIGURE 3M. - VARIATION OF PROPELLER DESIGN PARAMETERS WITH BLADE 
RADIUS FOR THE UNSWEPT SR-2 PROPELLER. 
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